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ABSTRACT 

Organic material found in meteorites and interplanetary dust panicles is enriched in D and 1:> N. This is consistent 
with the idea that the functional groups carrying these isotopic anomalies, nitriles and amines, were formed by 
ion-molecule chemistry in the protosolar nebula. Theoretical models of interstellar fractionation at low temperatures 
predict large enrichments in both D and 1: *N and can account for the largest isotopic enrichments measured in 
carbonaceous meteorites. However, more recent measurements have shown that, in some primitive samples, a large 
i5 N enrichment does not correlate with one in D, and that some D-enriched primitive material displays little, if 
any, l5 N enrichment. By considering the spin-state dependence in ion-molecule reactions involving the ortho and 
para forms of H?, we show that ammonia and related molecules can exhibit such a wide range of fractionation for 
both ,;> N and D in dense cloud cores. We also show that while the nitriles, HCN and HNC, contain the greatest 15 N 
enrichment, this is not expected to correlate with extreme D enrichment. These calculations therefore support the 
view that solar system 3:> N and D isotopic anomalies have an interstellar heritage. We also compare our results to 
existing astronomical observations and briefly discuss future tests of this model. 
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1. INTRODUCTION 

The isotopic enrichments measured in the primitive organic 
matter found in comets, interplanetary dust particles and me- 
teorites (Re mu sat et al. 2006; Floss el al. 2006; Aleon 2010), 
probably had their chemical origin in a low- temperature envi- 
ronment that facilitated ion-molecule isotopic exchange reac- 
tions, such as the dense interstellar medium (ISM) or the outer 
protosolar nebula (Mumma & Charnley 201 1: Messenger el al. 
2003). 

The enhanced D/H ratios seen in some primitive material 
(Messenger 2000; Robert 2003; Messenger et al. 2003; Keller 
ct al. 2004; Hartogh et al. 2011) almost certainly originated 
in low- temperature environments. Deuterium enhancement is 
initiated by the exothermic ion-molecule exchange reaction 
(Millar et al. 1989) going to the right 

H+ + HD -i H 2 D + + H 2 + 170K. (1) 

In dense cores where CO and other heavy molecules have 
condensed as ices on to dust grains (i.e,, become depleted), 
H 2 D + molecules can continue to react with HD, producing HD} 
and Dy (Roberts et al. 2003). This produces very high levels 
of deuteralion in the remaining gaseous species as well as, 
through the resulting high atomic D/H ratios, in molecular ice 
mantles (e.g., Ceccarelli et al, 2007). A similar D fractionation 
chemistry can also occur in the cold regions of proloplanetary 
disks (Willacy 2007). 

In both J upiter Family and Oort Cloud comets, U N/ 1:) N ratios 
are measured (~ 1 30—170) for both CN and HCN (Hutsemekers 
et al. 2008; Boekelee-Morvan el al. 2008) that are significantly 
lower than both the terrestrial and protosolar values (272 and 
440, respectively; Marty 2012). Primitive materials can exhibit 
large bulk 1:> N enhancements ( n N/ I:, N ~ 180-300), and high 
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spatial resolution measurements further show that they also 
contain distinct “hotspots’' which exhibit the greatest D and 
1:> N enhancements known ( 14 N/ J5 N < 100) (Messenger 2000; 
Aleon et al. 2003; Floss et al. 2006; Busemann et al, 2006; 
Briani et al. 2009a; Aleon 2010), 

Theoretical models of nitrogen isotopic chemistry in molecu- 
lar cloud cores, where gaseous molecules are being depleted by 
accretion on to dust grains, can reproduce the range of n N/ 1:) N 
ratios found in primitive matter (—50-280; Aleon 20 1 0) through 
exchange reactions between l:> N atoms and molecular ions such 
as i - | N 2 H + and HC I4 NH + (Charnley & Rodgers 2002; Rodgers 
& Charnley 2008a, 2008b). These models predict that N 2 , NH 3t 
HCN, CN, and HNC molecules should possess the greatest 
,5 N enrichments. There arc two ,5 N-fraciionating pathways 
(Rodgers & Charnley 2008b): a slow one to N 2 and ammo- 
nia MO 6 years) and more rapid one to HCN and other nitriles 
(^ ] (P years). 

Recent, albeit sparse, astronomical observations of starless 
cores indicate that interstellar NH 2 and N 2 are not as enriched 
as predicted for the ISM or as measured in the solar system: 
n NH 3 / I5 NH 3 * 334 (Lis etal. 2010) and m N 2 H7 h N I 5 NH + % 
446 (Bizzocchi et al. 2010). On the other hand, preliminary 
measurements indicate that interstellar HCN and HNC are 
highly enriched in ,5 N (Milam & Charnley 2012). These 
discrepancies arc potentially related to different fractionating 
pathways and timescales. 

However, the most pressing problem lor a direct ISM -solar 
system isotopic connection concerns the fact that, while ,:> N 
and D meleoritic hotspots do seem to correlate in some samples 
(e.g., Aleon el al. 2003), they clearly do not in others (e.g., 
Busemann et al, 2006; Courier et al. 2008). Yet, as has been 
noted by several authors (Alexander et al. 2008; Briani et al. 
2009b; Bonal et al. 2010; Aleon 2010), the low-temperature 
interstellar environments most conducive to producing large 
,:> N enrichments, should also produce concomitantly enormous 
molecular D/H ratios (e.g,, in ammonia or HCN). Thus, one 
would expect 1:> N and D hotspots to always be spatially 
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cor related, contrary to what is seen. This presents a serious 
challenge for ion-molecule fractionation chemistry. 

In cold interstellar environments, nuclear spin stales can be 
an important factor in determining chemical reaction rates (e.g.. 
Flower ct al. 2006). The major interstellar collision partner, 
molecular hydrogen, has a 170.5 K difference in zero-point 
energy between the para (anti parallel spins) and the ortho 
(aligned nuclear spins) forms; this allows some endoergic 
reactions to proceed, even at low temperatures, if a substantial 
fraction of the H 2 is in the higher energy form, o- H 2 . Pagan i 
et al. (201 1) demonstrated that a high H 2 ortho-to-para ratio 
(OPR) therefore inhibits (‘‘poisons”) the production of H 2 D + in 
reaction (1), suppressing the overall deuteration. 

The fractionation of b N may also be influenced by the OPR 
of H 2 . Ammonia formation is initiated by the production of N + 
from N 2 by He + , which then reacts in 

N + + H 2 — > NH + + H, (2) 

followed by a sequence of ion-molccule reactions with H 2 
through to NH4 and a final electron dissociative recombination. 
Reaction (2) has an activation energy of £200 K (Gerlich 1993) 
which, at low temperatures, can be overcome by the internal 
energy of 0 - H 2 (Le Bourlot 1991). Thus, ammonia formation 
in dark clouds is dependent on the H 2 OPR, a distinction not 
made in earlier fractionation models (e.g., Rodgers & Charnley 
2008b). A recent re-assessment of the OPR dependence in the 
original experimental data by Dislaire et al. (2012) indicates 
that previous work has overestimated the low-temperature rate 
coefficient for reaction (2), involving o- H 2 , by almost three 
orders of magnitude (cf. Le Bourlot 1991). Thus, the precise 
abundance of o- H 2 could play a pivotal role in producing a 
diverse range of D- b N fractionation in interstellar precursor 
molecules. 

In this Letter, we quantify the effect the H 2 OPR has on 
interstellar b N fractionation and demonstrate that this can 
account for the isotopic anomalies observed in primitive solar 
system materials. 

2. MODEL 

We consider the chemical evolution of the central regions of 
a static prestcllar core (Charnley & Rodgers 2002) with a gas 
density n( H 2 ) = I0 6 cm” 3 , a temperature of 10 K and a visual 
extinction of A\ > 10 mag. Cosmic-ray ionization occurs at a 
rate off = 3 x 10" 17 s"’ 1 . The chemical model used here is based 
on Rodgers & Charnley (2008a, 2008b) 2 and Figure I shows the 
most important reactions for b N fractionation. Helium, carbon, 
oxygen, and nitrogen are included at elemental abundances 
given by Savage & Scmbaeh (1996). We are interested in cores 
that arc just about to form protostars and so we assume that all 
carbon is initially bound up in CO, with the remaining oxygen 
in atomic form, and the nitrogen is partly atomic with 50% 
in molecular and atomic forms, respectively. 3 The elemental 
l4 N/ b N ratio is assumed to be 440, with the same fraction of 
b N in atomic form initially, b N/ b N l4 N = 1. 

The gas-grain interaction is modeled according to Charnley 
(1997); all neutral gas-phase species, except H 2 , He, N, and N 2 , 
stick and frcczcout onto grains upon collision. No grain surface 
chemistry or desorption is considered, except for slicking of 

2 We do not consider ilie neutral-neutral isotope exchange reactions 
proposed l>v Rodgers & Charnley (2008b). 

2 The effect of different N/N 2 ratios on this chemistry is discussed in detail 
by Rodgers & Charnley (2008a), 
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Figure 1. Chemical network showing the main reactions responsible for *'N 
enhance mem in nitriles and ammonia. 

H atoms (with an efficiency of 0.6) followed by reaction to 
produce H 2 molecules which are ejected upon formation (see 
Section 2. 1 ). 

2. /. Inclusion of Spin-state Reactions 

The reaction network has been expanded to contain ortho 
and para forms of H 2 . HJ, and H}, and couples 279 gas- phase 
species by 4420 reactions. Both spin types have been considered 
equally likely to react with a given species, with the exception of 
reaction (2), where new rates of Dislaire et al. (2012) for ortho 
and para H 2 , k 2 .„ = 4.20x 10 ” lo ( r / 300)" OJ7 exp(-44.5/T) 
cm 3 s“ J and £ 2 ./> = 8.35 x 1()“ 10 exp(— 168.5/T) cm 3 s'" 1 , 
respectively, have been adopted. Other temperature-dependent 
H 2 reactions have energy barriers that are too high lor the 
ground -state energy difference between the spin types to have a 
significant effect at 10 K. 

Initially, the H 2 OPR is set to the high- temperature statistical 
equilibrium value of 3. Measurements in the diffuse ISM by 
Crabtree et al. (201 1) are consistent with an H 2 OPR around 1 ; 
wc find adopting this lower value has little effect on our results. 
The reaction 

H + + o- H 2 — > H + + p~ H 2 (3) 

is exothermic by 170.5 K and so o-H 2 is efficiently convened 
into p- H 2 at low temperatures. We assume that the H 2 formed 
and ejected from grains has an OPR of 3; this ensures a 
steady- si ale OPR is reached. Reaction rates for conversion 
between ortho and para species, as well as for their interaction 
with electrons and cosmic rays, have been implemented from 
Walmsley et al. (2004). Branching ratios for o/p-Hy formation 
from Hj and H 2 , and rates for Hy interaction with CO and N 2 , 
are taken from Pagan i et a). (2009). Furthermore, all reactions 
involving proton transfer from H* to neutral molecules have 
been assigned branching ratios according to Oka (2004) and no 
spin conversion is assumed to lake place in reactions involving 
electron transfer to H| . Apart from these special cases, reactions 
arc assumed to produce ortho and para versions of H 2 , H£ and 
Ht with equal probability. 
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Figure 2, Left panel: lime evolution of ihe nitrogen chemistry in dense cores, compared 10 CO and o- lb. Crucial limes for ihc ,3 N fraeuonaiion in nitrogen hydrides 
are marked by vertical dolled lines, labeled A, B, and C. Right panel: evolution of the u N/ 1:, N ratio in gas-phase molecules and bulk of ammonia and nitrile ices. Black 
curves represent the current model, while red curves show the corresponding isotopic ratios in the mode! of Rodgers & Chamley (2008a). The assumed elemental value 
of 440 is marked by a dotted horizontal line. Values for gas- phase nitriles are omitted after x 10 s years when abundances for both isoiopologues are essentially 
zero due to freezeout. Dashed curves represent ice- phase species. 


2.2. Other Reactions 

Since nitrogen fractionation is initiated by He + reactions 
(see Figure 1), the He + chemistry dictates the timescale of 
fractionation. In particular, the rate for the reaction 

He + + H 2 — > H + + H + He (4) 

has previously been assigned widely different values at 10 K, 
ranging from 1J x 10 ' 15 cm 3 s ' 1 (Rodgers & Charn ley 2008b) 
to 2.5 X 10 13 cm 3 s" 1 (Walmsley et al. 2004). We have adopted 
the measured low-temperature rate k+ — 3.0 x 10 _N enr s 1 
(Schaueretal. 1989). 

In order to qualitatively compare the temporal patterns in l;> N 
fractionation to the expected levels of deute ration during core 
evolution, we also include the chemistry for ortho and para 
H 2 D + in the model since o- H 2 can overcome the barrier for the 
reverse of reaction (1). i.e.. 

H 2 D + + oH 2 — ► Ht + HD. (5) 

The most important reactions for spin- stale -dependent o/p - 
H 2 D + formation and destruction tire adapted from Walmsley 
el al. (2004), and we assume a deuterium abundance corre- 
sponding to HD = 1.5 x 10“^ relative to elemental hydrogen. 

3. RESULTS 

Figure 2 shows the evolution of the major chemical species. 
The ammonia abundance builds up to about 10“ 7 with respect 
to H 2 around a core age of 10 5 years (left panel), which is 
comparable to what is typically observed (Hol/.el el al. 2004). 
Similarly, the peak HCN and HNC agree well with values 
observed in dense cores (Padovani et al. 201 1). At this time 
the CO is depleted by a factor of > 1()\ A steady-state H 2 OPR 


of 2.3 x 10' \ consistent with attempts to quantify the OPR 
in the cold dark cloud B68 (Troscompt et al. 2009), is reached 
after about 10 6 years. 

The right panel of Figure 2 shows the evolution of the 
molecular n N/ 15 N ratios in the current model, compared to 
the corresponding results from the Rodgers & Chamley (2008a) 
model. As expected (see, e.g., Figure 1), the inclusion of a 
spin-state dependence does not significantly affect the evolution 
of fractionation in the nitriles: the i:> N enhancement starts 
to buildup from the initial pool of atomic 1<S N early on, and 
continues to increase until gas-phase abundances have dropped 
to below 10” 13 . The ]4 N/ ,5 N ratios in gas-phase HCN and HNC 
lie in the range ~50-260. When all nitriles are frozen out, the 
Iy N enhancement in the bulk of nitrile ices (CN, HCN, HNC, 
and H 2 CN) is —5. 

The inclusion of spin-slate dependence does however have a 
dramatic effect on (he 1:> N fractionation in ammonia. The relative 
delay in the timing of the late enhancement peak in the current 
model, as compared to Rodgers & Chamley (2008a), stems 
from the higher rate 0-3 x) of reaction (4) (see Section 2.2). 
However, the main difference is that the NH3 enrichment is 
not monotonic in time. The sharp decrease starting around 2 x 
1CP years (time B) (i.e., an increase in J4 NH;*/ ,;> NH;0, is caused 
by the inefficiency of reaction (2) as the c?-H 2 abundance drops 
(see Figure 2). 

To understand what causes the fractionation patterns, we 
compare abundances to i so to po log ue ratios at crucial times in 
the evolution, marked by vertical lines A. B, and C in Figure 2. 
Initially, the limiting reaction for fractionation is the production 
of J: 'N + , and the l4 N/ l:> N ratio in ammonia is simply double 
that of n N 2 / 14 N 1, 'N. As gas-phase CO gets depleted after 
'"HO 4 years (A), more He + becomes available to produce b N\ 
which starts to enrich the nitrogen hydrides beyond that of N 2 . 
However, after — KB years, the H 2 OPR drops below 2 x 10“ 7 
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Figure 3. Nitrogen fractionation in ammonia and HCN, both gas-phase (solid) and bulk ice (dashed) in blue curves according io left-hand axis, compared lo deuterium 
fractionation in H 2 D + /Ht (green curve and scale on right-hand axis). enhancements are given relative to the elemental l? N/ l l N ratio of 1/440. while D 
enhancements are relative to elemental D/H = 1.5 x 10"*. 


(13), reaction (2) becomes less efficient, and more ,5 N + is instead 
circulated back into molecular nitrogen by 

l5 N + + !4 N 2 r=i m N + + 14 N J "N (6) 

(Terzieva & Herbsl 2000). lowering the b N fraction in ammonia 
to less than half the elemental fraction. The reason for ammonia 
to subsequently become ,:) N enriched after 10 6 years (C) is 
that the molecular nitrogen fraction decrease substantially. This 
both reduces the fraction of l5 N + converted back into molecular 
form (Figure 1) and decreases the amount of 14 N + available to 
produce l4 NH.v 

Fractionation changes in the bulk of amine ices follows the 
pattern of gas-phase fractionation, but with a time lag and not 
as pronounced extreme values. While the bulk of amine ices on 
grains has ,:> N fractions within 30% of the elemental fraction, 
at any given time the fractionation in the outer monolayer may 
show much larger variations (Rodgers & Charnley 2008a). 

In Figure 3 the calculated deuterium fractionation in H| 
is plotted together with the molecular ,:> N fractionation. The 
D/H ratios in other gas-phase molecules, in atoms, and therefor e 
also in the accreted molecular ices, will closely follow that of 
U*. The increase in deuteration follows very closely the drop 
in abundance, and after 2 x 10 5 years (time B), ratios of 
D/H > 1% are predicted. These trends in both deuteration and 
OPR evolution agree well with previous investigations (Pagani 
et al. 20 1 l ; Flower el al. 2006). 

Figure 3 shows that gas-phase chemistry can indeed repro- 
duce the requisite range of D and enrichments in nitrile and 
amine functional groups. Three fractionation patterns develop 
as the core evolves as following. 

I. For times less than about ^2 x HP years, gaseous nitriles 
can develop large enrichments in ln N (HC h N/HC ,:> N ^ 
50-260) but will only be modestly enriched in D, as- 
observed (e.g., DNC/HNC ^ 0.01-0.05; Hirota el al. 
2001); this will be reflected in the accreted ices. Ammonia 
will be similarly enriched in D but less enriched in l:> N than 


the nitriles. At this stage the nitrile gas-grain fraciionation 
is complete. 

2. Between x I0 3 -10 6 years, the remaining ammonia can 
become highly fractionated in D but is either not enriched, 
or even depleted, in b N. 

3. Beyond 10 6 years, ammonia becomes very highly enriched 
in both D and l5 N. 

Thus, the model naturally produces a range of D- ] *N frac- 
tionation in nitrile and amine functional groups that is quali- 
tatively consistent with that found in primitive organic matter. 
In general, the highest l:> N enrichments and lowest D enrich- 
ments should be carried by nitriles. Amines could account for 
the fractionation of material that is highly enriched in D only, 
as well as the material that is simultaneously highly enriched in 
both D and l5 N. This functional-group dependency has already 
been noted in analyses of meteoritical samples. However, an im- 
portant caveat is that conversion of nitriles to amines by H and 
D additions on grain surfaces (c.g., hydrogenation of HCN to 
CH 2 NH and CH 3 NH 2 ; Theule et al. 2011) has to be negligible. 
The 30% enrichment of the bulk ammonia ice cannot account 
for all the enhanced ammonia fractionation measured in car- 
bonaceous chondrites ( 14 NH 3 / I: *NH;, % 185-256; Pi/zarcllo & 
Williams 2012). However, we predict I4 NH 3 / IS NH 3 % 100 as 
the highest gas-phase enrichment and this will lead to higher 
,4 NH 3 / ,5 NH 3 ratios in the uppermost monolayers of the strati- 
fied ices (Rodgers & Charnley 2008a). 

Interstellar gas-phase l4 N/ I5 N ratios in nitriles have been 
measured in several Galactic sources selected for being unlikely 
lo have experienced significant isotopic fractionation (Adande 
& Ziurys 2012). For cold dark cores, fewer measurements 
exist; our model predicts very low n N/ 15 N ratios for HCN and 
HNC. For ammonia, the observed J4 NH 3 / ,5 NH 3 ratios in dense 
protostellar cores (^340; Lis et al. 2010) are consistent with 
the modest enrichments which can occur within lO- years. The 
measurement of n NH 3 / I:> NH 3 > 700 in L1544 (Gcrin et al. 
2009) is consistent with the prediction that b N nuclei can be 
underabundant in ammonia. It also suggests a chemical age of 
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^3 x I0 5 years for the fractionated ammonia gas (Figure 3), 
and that it should exist in a region of the Li 544 core spatially 
distinct from HCN and HNC (cf. Padovani et al. 201 1). 

4. CONCLUSIONS 

By considering the ortho-para dependence in ion-molecule 
reactions involving H 2 , ammonia, and related molecules are 
expected to exhibit a wide range of fractionation in both ,;> N and 
D, consistent with those found in primitive solar system organic 
matter. These include very large enrichments in both D and 
]5 N, as found in previous models, but also low l5 N enrichment 
or even depletion, coupled with modest to large D/H values. 
Nitrile functional groups on die other hand, while being the most 
likely source of b N hotspots, are not predicted to correlate with 
extreme D enrichment. These calculations support (he premise 
that the ,5 N isotopic anomalies found in meteorites were set in 
functional groups formed in a cold molecular cloud. 

The predicted ,4 N/ I:> N and D/H ratios can account for the 
limited measurements existing for interstellar cloud cores. To 
further test and constrain these theoretical models, surveys for 
key molecules and measurements of relevant 14 N/ 1; "N ratios in 
molecular clouds are required. A more extensive analysis of how 
correlations between deuterium and ,3 N enrichment in specific 
molecules are affected by spin -state-dependent chemistry is be- 
yond the scope of this Letter and will be considered elsewhere 
(E. S. Wirsirom et al. 2012, in preparation). Predictions from 
these models can be tested by measurements of various mul- 
tiply substituted isolopologues (c.g., DC^N, 1:, NH 2 D; Gerin 
et al. 2009), which will be facilitated with the Atacama Large 
Millimeter Array. 
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